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ly-j ■ We present a systematical analysis of uncertainties in the helioseismological 

determination of quantities characterizing the solar structure. We discuss 
the effect of errors on the measured frequencies, the residual solar model 
dependence and the uncertainties of the inversion method. We find Y p h = 
0.238 - 0.259, R b /R & = 0.708 - 0.714 and p b = (0.185 - 0.199) gr/cm 3 (the 

£NJ ■ index b refers to the bottom of the convective envelope). In the interval 

0.2 < R/R & < 0.65, the quantity U=P/p is determined with and accuracy 

Q\ . of ±5°/ 00 or better. The predictions of a few recent solar model calculations 

\ are compared with helioseismological results. 

Of 

S ■ I. INTRODUCTION 

+- > ■ 

A Standard Solar Model (SSM) may be defined as a description of the solar interior which 
reproduces the observed properties of the sun within observational errors, by adopting a set 
of physical and chemical inputs chosen within the range of their uncertainties. 

From the point of view of stellar evolution, the sun is a well known structure. Confidence 
in SSMs rests on the success of stellar evolution theory to describe many, and more com- 
plex, evolutionary phases in good agreement with observational data. Before the advent of 
helioseismology, however, a solar model had three (essentially) free parameters - the initial 
helium and heavier elements abundances, Y and Z, and the mixing length parameter a - 
for producing three observables: the present radius, luminosity and heavy elements con- 
tent of the photosphere. From this point of view, the success of SSMs to reproduce these 
observables may look as a not too big accomplishment. 

The solar neutrino puzzle (see e.g. |H§), i.e. the disagreement between the SSM pre- 
dictions and the results of the four solar neutrino experiments 0-0) originated many dis- 
cussions about the validity of the SSMs and prompted many attempts to build alternative 
solar models, none however capable of solving the discrepancy. 

In recent years helioseismology has added important information on the solar structure. 
A SSM has now to account for these additional data. 
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Helioseismology (see || for an introduction and flJTOl for extensive reviews) provides 
severe constraints and tests of solar model calculations. For instance, from inversion (decon- 
volution) of helioseismological data one can infer the photospheric helium abundance Y p h 



|TT| , p!2[| and the depth of the convective zone Rb [y|. The helium mass fraction, which was 
used as an essentially free parameter to match the observed solar luminosity, is now strongly 
constrained. Still one can adjust a in order to get the proper solar radius, R®; however there 
is no parameter which can be tuned to get Rb- The comparison between the theoretically 
predicted value of Rb and the value inferrred from helioseismology is a test of the model. 

Helioseismology can tell a lot about the solar interior, however, it cannot replace SSMs. 
As an example, helioseismology determines the sound speed down to the solar core to high 
accuracy, but it cannot give the temperature profile in the energy production region. In 
order to calculate the several neutrino luminosities (pp, 7 Be, CNO, 8 B...), which crucially 
depend on this profile, SSMs are needed. 

For people interested in solar neutrinos, the basic question is thus: 

• Which solar models - if any - are acceptable, i.e. consistent with helioseismology? 

Before addressing this question, one has to provide a quantitative answer to the following 
one: 

• How accurate are the helioseimological determinations of quantities characterizing the 
solar interior? 

The present paper is essentially devoted to this issue: we perform a systematic and 
possibly exhaustive investigation of the uncertainties of the helioseismological approach, in 
order to estimate the global error to be assigned to helioseismological determinations of the 
solar properties. With this spirit, we analyse the helioseismological determinations of several 
physical quantities Q characterizing the solar structure. Concerning the outer part of the 
sun, we discuss Y ph , R b , the sound speed q, and density pb at the bottom of the convective 
zone. Then we consider the "intermediate" solar interior (x=R/R® = 0.2 — 0.65), analysing 
the behaviour of the squared isothermal sound speed, U —P/ p. Finally we investigate the 
inner region (x < 0.2), where nuclear energy and neutrinos are produced. 

For each quantity Q we determine the partial errors corresponding to each uncertainties 
of the helioseismological method. This approach - which will be better elucidated in sections 
1 1| and |H| - is needed since helioseismology measures only the frequencies {Q} of solar p- 



modes, and quantities characterizing the solar structure are indirectly inferred from the 
{f2}'s, through an inversion method. Schematically, the procedure is the following: 

a) One starts with a solar model, giving values Q mo d and predicting a set {fl mo d} of 
frequencies. These will be somehow different from the measured frequencies, Q® ± AQ® 

b) One then searches for the corrections q to the solar model which are needed in order 
to match the corresponding frequencies {Vt mod + u} with the observed frequencies {O }. 
Expression for u = w(q) are derived by using perturbation theory, where the starting model 
is used as a zero-th order approximation. The corrections factor q are then computed, 
assuming some regularity properties, so that the problem is mathematically well defined 
and/or unphysical solutions are avoided. 
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c)The "helioseismological value" Q® is thus determined by adding the starting value and 
the correction: 



Qe = Qmod + q- (1) 

There are three independent sources of uncertainties in this process: 

i) Errors on the measured frequencies, which - for a given inversion procedure - propagate 
on the value of Q Q . 

ii) Residual dependence on the starting model: the resulting Q Q is slightly different if one 
starts with different solar models. This introduces an additional uncertainty, which can be 
evaluated by comparing the results of several calculations. 

iii) Uncertainty in the regularization procedure. Essentially this is a problem of extrap- 
olation/parametrization. Different methods, equally acceptable in principle, yield (slightly) 
different values of Q Q . 

It has to be remarked that, in view of the extreme precision of the measured frequencies, 
AQq/Qq ;$ 10~ 4 [I^-|I71, uncertainties corresponding to ii) and iii) are extremely important. 



We consider all these uncertainties, as well as the possible ways of combining them to 
estimate the global error. 

All this will provide a framework for a quantitative comparison between solar models 
and helioseismology. In the final part of the paper we test the predictions of a few recent 
(standard) solar model calculations. We feel that similar test should be performed, when 
alternative solar models are proposed. 

The paper is organized as follows: 

A short review of the inversion methods is given in section [IT], where we also define the 
free parameters and determine their allowed ranges. The individual uncertainties mentioned 



above are operationally defined in section |HI| , where possible definitions of the global errors 



are also discussed. In section |I^ the properties of the outer sun (photospheric helium 
abundance and quantities characterizing the bottom of the convective zone) are examined, 
and results are compared with previous analyses. Section [V| is devoted to the intermediate 
region (x = 0.2 — 0.65) whereas section [VT] concerns the energy production region. In section 
|V11| the results of a few recent solar models are compared with the available helioseismological 
constraints/tests. 

The concluding section describes the main results of the paper, summarized in Figs. |2| 
and ^| and in Table 0. 



II. HELIOSEISMIC DATA AND INVERSION METHODS: A SUMMARY 

a) Observational data 

The helioseismic data that have been used for probing the solar interior structure are 
frequencies of acoustic oscillations (p-modes) corresponding to Y® spherical harmonics, with 
the degrees / ranging from to some 150. At low degrees there are typically 15 modes 
corresponding to rather high radial degrees, n > 10. The total number of modes which 
has been used in structural inversions exceeds two thousands, which is only a small fraction 
of the total number of modes with measured frequencies. The latter number is higher by 
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two orders of magnitude. However, modes corresponding to I > 150 are not really useful 
for probing the structure of the sun's interior, because they propagate only in the shallow 
outermost layers and also because the accuracy of the corresponding frequencies is much 
worse. Tesseral (m 7^ 0) mode frequencies yield an information about the sun's internal 
rotation and magnetic field. In the present application, the only relevance of these data is 
to justify ignoring effects of the centrifugal and Lorentz forces on the radial structure. 

The most exploited helioseismic data set is one from measurements at the Big Bear Solar 



Observatory in 1986 (BBSO86 This set contains no data for modes with I < 3 and has 

been supplemented with data on I = 0-3 modes from various sources. The most extensive 
data on such modes is from the network of automatic telescopes operated by the group of 
Birmingham University (BISON [|15|]). Recently, data covering the whole range of /-values 
became available from the LOWL instrument |16[ and from the GONG network |TJ 

b) Expression for w = u>(q) 



The basic assumptions made in seismic probings of the sun's interior structure include 
hydrostatic equilibrium and adiabaticity of oscillations. In early inversions an asymptotic 
approximation for the radial eigenfunctions was used leading to a simple integral equation 
connecting the adiabatic sound speed, c(R), to oscillation frequencies f^ jTl . However, this 
approximation is inadequate, especially for probing the deep interior and in most of the 
subsequent works numerical solutions have been adopted. The price is a need to use a 
SSM as a starting model about which the hydrostatic equations are linearized. There is 
an implicit assumption that the model provides a sufficently close description of the sun's 
structure. With this assumption, the variational principle following from the equation for the 
adiabatic oscillations, leads to an integral equation connecting the differences between the 
solar and the model functions describing radial structure (e.g. Ref . |18| ) to the corresponding 
differences in mode frequencies. The equation, however, cannot be directly applied to infer 
the difference in the structural functions because the assumption of hydrostatic equilibrium 
of the mean structure as well as that of adiabaticity of oscillations break down near the 
surface. To account for these departures an ad hoc term of the form F(Q)/I, where Q is the 
frequency an / is a suitable moment of inertia, must be added . Such a form follows from 
the fact that there the radial eigenfunctions are nearly /-indepedent. The resulting formula 
may be written in the following form, 

» ^£"mod ' j J tymod •> mod *j 

where j = (l,n) identifies the mode, Q(x) is a structural function, T= (cflogP /dlogp) a d is 
the adiabatic exponent, x = R/R & and: 

q = Qq — Qmod (3) 

7 = — T mo d 

There is freedom in choosing Q. It could be density p, pressure P, or any combination 
of these quantities and their derivatives. All such functions are connected through the 
linearized mechanical equilibrium condition. 
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Even with a very large number of Eqs. @, equal to the number of all modes in the data 
set, a determination of the three functions: q(x), j(x) and F(Q) is not possible without the 
additional assumption that the unknown functions vary slowly. One may eliminate one of 
the functions by making use of a thermodynamical relation T = T(P, p,Y), where Y is the 
fractional He abundance, and of two, well justified simplifications, which are: 

(1) Y=Y p h in the convective envelope, and 

(2) 7=0 in the radiative interior. 

With these additional constraints the unknown function 7 is related with the unknown 
number Y® h . This approach was first adopted in [IT5| with the choice that the single directly 
determined function was U = P/p. 

Thus Eqs. @ are written in terms of the unknown function 

u(x) = U® - U mod , (4) 

the function F(Q) and of the unknown number y p h = Y® ft — Y™fi d . One advantage of this 
formulation of the inverse problem is that one obtains directly the surface abundance of 
He. We follow this approach here. Inversion without reference to thermodynamics has been 
done by Antia and Basu . 



c) The regularized least square method and its free parameters 

One way of making use of the assumption about the slow variability of u/U mo d and 
F(Q) is a discretization in terms of known functions. Following |TJJ, we use cubic splines 
in the first case and Legendre polynomials in the second case. The coefficients in these two 
representations together with y p h are determined by the least-squares method. A simple- 
minded application of this method results in a solution for u/U mo d which exhibits artificial 
oscillations. A cure to this problem is regularization (see e.g. |2~1|| ) which consists in adding 
to the usual \ 2 an additional term. Here we minimize 

(5) 

where J is the total number of the frequency data, AQ Q are the errors on the measured 
frequencies and A is a control parameter. This form of regularization is referred to as the 
second derivative smoothing. Alternative forms of regularization and influence of the choice 
of A were investigated in . 




The regularization parameter A was chosen to be the minimum value that still suppresses 
the oscillations. Numerical experiments were conducted in which attempt was made to 
reproduce known differences in u(x) between two solar models from differences in p-mode 
frequencies. The set of modes was the same as in the observational data with weights 
determined by the observational errors. The conclusion from these experiments was that for 
x > 0.1 it is possible to reproduce u(x) very accurately and independently of the choice of 
the type of the regularization. The results are also insensitive to A in a wide range of values. 
However in the inner core (x < 0.1) a regularization needed to avoid the oscillatory behaviour 
always flattened the real steep increase in u/U mo d- It became clear that the method is not 
applicable for probing this part of the sun. The reason is that only a small fraction of modes 
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in the data sets exibits any measurable sensitivity to the sound speed changes in the inner 
core and the effect of regularization dominates over the data. 

Uncertainities on j p h and u(x), as well as on other indirectly determined structure func- 
tions, following from observational errors, are best estimated by random number simulation 
having normal distribution with the half- width given by AQ & . There is an uncertainty 
following from the freedom in: choice of the regularization parameter A, number of bins, 
rib, and number of the Legendre polynomials, rip. This has been discussed in |T9LE2LE3 



Changes in rib m ay always be compensated by changes in A and will not be discussed here. 
Admissible ranges for A and rip follow from considering the behavior of Xreg- An increase in 
A, as seen in Eq. (|j), causes an increase in xl eg y the same effect being caused by the decrease 
of tif- Specifying the maximum departure from the minimum value XrJg MI we can se ^ an 
upper bound for A and a lower bound for np. Concerning A, one adds the requirement of 
avoiding artificial oscillations. Dealing with the real data we may still distinguish the arti- 
ficial oscillations from features implied by data. The former have the half-wavelength equal 
to the bin length and their amplitude increases sharply with the decrease of A. Appearence 
of the oscilllation sets a lower bound for A. 

Let us summarize the role of the parameters of this method, which will be referred to as 
the regularized least square (RLS) method: 

i) There are two parameters, A and rip. 

ii) The seismic model of the sun is obtained by taking A=0.001, as this is the minimal value 
for which unphysical oscillations are avoided, and rip=20 as this corresponds to the minimal 
degree of the polynomial such that xleg * s stable. 

iii) For this set of parameters, by using the "model S" of Ref. |24j (hearafter referred to as 
JCD), as a starting model, we obtain xle g MIN ^ = ^-^i which provides an indication of the fit 
quality. 

iv) For the parameters A and rip we consider as acceptable those values which yield a xfeg 
not substantially degraded with respect to the minimum. As a prescription borrowed from 
statistics we accept the parameters if Xreg — Xreg IN ^ + 1- Numerically we found that this 
implies A < 0.1 and riF>10. 



d) Methods and parameters for the solar core 



As already remarked, the method described above cannot be used for probing the core 
within x—0.1. However, a fairly accurate seismic sounding of that region should be possible. 

It was found [23] that frequencies for some 30 solar p-modes were changed by more 
than AQ & when U in this innermost part of the core was changed by just one percent. 
For that region the optimal localized averaging (OLA) method is a better choice. In this 
method we do not try to determine functions. Instead, we try to determine mean values 
weighted with (possibly narrow) kernels centered at selected {a;o} values. We consider a 
linear superposition of individual kernels for the modes present in the data set, 

j 

K u (x ,x) = Y, c j( x o)£ u> j{x) (6) 

3 

and determine the coefficients Cj at selected points. An application of the classical method 



6 



of Backus & Gilbert [f25| to helioseismic data was described by Gough & Thompson |L8 . 



The subtractive optimal localized averaging (SOLA) method adopted in [g^l an d used also in 
this paper was developed by Pijpers & Thompson [PEfl . Here one tries to construct kernels 
K u (x ,x) which are as close as possible to Gaussians, G, centered in x Q and characterized 
by their half- widths at half- maximum, w. Subject to minimization is the quantity 



K u — G 



X — Xq 



w 



2 J 



dx + nJ2^ n lr ( 7 ) 

j 



The second term, with a trade-off parameter /i, is added to avoid large error magnification. 
For a specified w, larger [i leads to smaller errors in the localized mean but the kernels may 
differ significantly from Gaussians and therefore w cannot be regarded as the measure of 
localization. Once the coefficients Cj are determined the localized averages are: 
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A comparison of results of inversions with the two methods, given in p2| , |23| , shows a very 
good agreement everywhere except for the inner core. Results obtained with the SOLA 
method are more realistic because it does not introduces an artificial smoothing. Thus, 
for x > 0.1 we may rely on the functional form of u(x) as obtained by means of the RLS 
method. However, for a reliable and accurate probing of the inner core a different approach 



is required. In a hybrid method of seismic model construction developed in [22], the RLS 
method is used to determine F(Q), y ph and u/U mo d in the [xj, 1] interval, whereas the SOLA 
method is used to determine a single average value in the inner core around the point Xq. 
The effect of F(f2)and y p h is removed from the data. The function u/U mo d in the [0,Xf] 
range is then represented as a three term power series in x 2 thereby satisfying the boundary 
condition at x = 0. The three coefficients in the series are determined using the average 
SOLA value with the true kernel around xq and the continuity conditions at x — Xf for 
u/U mo d and its derivative. The method introduces its own uncertainties resulting from some 
freedom in choosing xq and Xf. Regarding the role of the parameters of the hybrid method 
we point out that: 

i) The constraint on Xf is that the point should be located in the region where the SOLA 
values agree within the observational errors with the values inferred by the RLS method. 
This occurs for Xfu > 0.075. We shall explore the effects of varying Xf in the range [0.075 
-0.125]. 

ii) We shall also study the effect of varying xq, while Xf is kept at the central value Xf=0.1. 
The range for x is suggested by two conditions: at too small values there is not enough 
information. On the other hand, when xq is too large and gets close to Xf there is no arm for 
extrapolating the information to small distances. As a consequence, somehow empirically, 
we chose xo in the range [0.04 - 0.06] and we adopt x = 0.05 as the "best" value. 

iii) The choice of the trade-off parameter, /i, for the hybrid method should be such as to 
minimize errors in u/U mo d implied by the measurement errors; we chose /i=0.001. 
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III. PARTIAL UNCERTAINTIES AND GLOBAL ERRORS 



We shall estimate how quantities Q characterizing the solar structure are changed when 
some input of the helioseismological determination are varied, one at a time. As a reference, 
we consider the values Q & obtained with the seismic model described in the previous section: 
frequencies are kept at the central measured values, the JCD model is used as the start- 
ing model, and the inversion parameters are kept at the "best" values A = 0.001, np=20, 
Xf it =0.1 and x =0.05. (The results of the seismic model are presented in Table [TI], second 
column, and Table [TV] again second column.) 

1) Observational errors . We consider the effect of varying the frequencies from the central 
measured values, within their observational errors. We took random choices of frequencies 
within ±Af2 and evaluated the corresponding values of Q. In this way we found maximal 
(minimal) acceptable values, Q sup (Qi n f)- We define the allowed range as: 

((5Q)exp Qsup Qinf (9) 

2) Starting model. For the central values of the observed frequencies and the best values 
of the inversion parameters, we study the differences resulting from using different starting 
models. We remind that JCD includes diffusion of helium and heavier elements, uses the 
Livermore equation of state and the Livermore opacity tables. To assess the relevance of 
these inputs, we performed inversion by using: 

a) a model (SUN241) with the same equation of state and opacities, but without diffusion 



( modified version of "model 0" from Ref. 23 



b) a model (SUN24) with the same opacities, the MHD equation of state and without dif- 
fusion ("model 0" from Ref. f23|] ). 

Exactly as before, for each quantity Q we evaluated the allowed range, (SQ) moc i- 



3) Inversion method . As explained in the previous section, two different inversion meth- 
ods are adopted, depending on the value of the radial coordinate. For x > 0.1, we use the 
RLS method and consider the effect of varying A and rip separately. The parameter A is 
varied in the range 0.001-0.1, see section [TJ. The resulting variation (maximum minus mini- 
mum in the explored range) induced on the quantity Q will be denoted as (5Q)\. The other 
parameter rip is varied in the range 10-20, see again section [H], and the resulting variation 
is denoted as (SQ) nF . 

For the inner region we use the hybrid method: the function U(x) is determined by 
the requirement that an interpolating polynomial matches the results of the SOLA method 
around x with that of the RLS method at Xf. As discussed in section ||, we consider Xf in 
the range [0.075 - 0.125] and xq in the range [0.04 - 0.06]. 

Although we are confident that our analysis is rather exhaustive, there is no well defined 
rule as to translate the differences (SQ)k into a global error AQ. Let us present two, somehow 
extreme, attitudes: 

a) The statistical approach. If one believes that, for each parameter, all possible values 
lie in the interval Q 1 ^ — Q^, one can interpret ±l/2(5Q)k, as a partial error. Furthermore, 
one can assume that the partial uncertainties somehow compensate each other, and use 
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the combination rule for the statistical fluctuations of independent variables, which (for 
Poisson-Gauss distribution) yields: 

(AQU = 4Jem ( 10 ) 

b)The conservative approach. May be that the parameter variation was not exhaustive, 
and what we found as extrema are not really so, but actually are quite acceptable values. In 
view of this, let us double the interval we found and interpret ±(5Q)k, as partial errors. Fur- 
thermore, let us be really conservative assuming that errors add up linearly. In conclusion, 
this gives: 

(AQ con ) = ±J2\(SQ)k\. (11) 

k 

In the following sections, we shall present both error estimates. The first one is rea- 
sonable when most of the error corresponds to the statistical uncertainty of the measured 
frequencies. Actually in our case the other errors (from the residual model dependence 
and/or inversion method) are generally dominant. Clearly they are of a more subtle nature, 
a sort of systematic errors. For this reason, we feel more confident with the conservative 
error estimate, Eq. flllD . 

Before concluding this section, let us briefly consider the question of the error correlation, 
which is important when considering the helioseismological implications on two quantities 
and (e.g. Y p h and R{). If the respective errors are AQ^ and AQ^ , we shall 
consider a solar model as acceptable if the predictions lie within the rectangle ±AQ^\ 
± AQ( 2 \ i.e we shall treat errors as uncorrelated, whereas actually there is a correlation: 
as an example, the extremum + AQ^ most likely corresponds to a set of parameters 
which does not allow - say - — AQ^ 2 \ The helioseismologically allowed area should be 
an ellipsis contained inside the rectangle . By choosing the rectangle as the allowed area, 
one is thus overestimating the error, which again corresponds to a conservative attitude. 



IV. PROPERTIES OF THE CONVECTIVE ENVELOPE 

A. The photospheric Helium abundance 

Knowledge of the helium abundance in the sun is of fundamental importance both to 
cosmology and to solar structure theory. Unfortunately, it is impossible to measure it, even 
in the atmosphere, by means of standard spectroscopy. A helioseismological method of 
determining Y p h was first suggested by Gough ]27| and developed in Refs. |28j and |p9| . The 



value of Y ph inferred for the sun, however, was not published, because it was feared that 
it might be too sensitive to systematic errors in the data that were available at the time. 
In \Y1\ an alternative method was developed. Unlike the former one, it does not rely on 



the validity of p-mode asymptotics. Instead, it requires the validity of linearization about a 
reference model. 

Several helioseismological determinations of Y p h have been published since 1991. The 
various estimates, which are collected in Table Q, often differ among each other by more than 
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the quoted errors, see Fig. |I[ Actually, these latter often reflect only the uncertainties of the 
measured frequencies, whereas the extracted value of Y p h depends on the inversion method 
and on the physical inputs, see the discussion in Refs. |3(J and Most of the estimates 
in Table | rely on the MHD equation of state. 

By omitting the highest value by Dappen et al. |32[ (the same authors later published 



significantly smaller values the total range of helioseismological determinations is: 



Y 



ph 



0.226 -0.260. 



(12) 



Our reference estimate is Y2 l =0.249 (we remark the we rely on the Livermore equation of 
state, which should be more accurate than MHD). Concerning the allowed range, the result 
of our conservative analysis, shown in Table [H], implies: 



Y, 



ph 



0.238 -0.259 



(13) 



which substantially overlaps with the previous range. From the same table one sees that 
experimental errors are of minor relevance, the value of Y p h being manly sensitive to the 
choice of the starting model; as discussed in Refs. |3C] and 
state is dominant. 



the effect of the equation of 



B. The bottom of the convective zone 



As emphasized in Ref. ||13||, the transition of the temperature gradient between being 



subadiabatic and adiabatic at the base of the solar convective zone gives rise to a clear 
signature in the sound speed. Helioseismic measurements of the sound speed therefore 
permit a determination of the location of the base of the convective zone. 
In Ref. |T3[ the following ranges were derived: 



R b /R G = 0.710 - 0.716 



(14) 



c b = (0.221 - 0.225) Mm/s, (15) 

where the quoted intervals include the uncertainty resulting from the inversion technique. 
The estimate of R b has been confirmed in Refs. [|33| , |30|1 . 

Actually, within the present uncertainty Rb and q, are not independent | I3]| . As well 



known, the lower part of the convective zone is very close to being adiabatically stratified, 
the adiabatic exponent T being close to 5/3, hence 

P oc p 5/3 . (16) 

Furthermore, throughout the same region, in the hydrostatic equation one can approxi- 
mate M x (the mass within x) with the total mass: 



dP _ GM Q 

dx R Q x 2 



P (17) 
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In terms of U = P/p the two equations give: 



dU _ 2 GM Q 

dx 5 Rq x 2 

and consequently 

2GM Q 



P (18) 



U{X) ~ 5"r7 (X_1) (19) 



which shows that U (hence c=y5/3£7) is uniquely known as a function of x. In view of 
the approximations which have been used [M x ~ M , T ~ 5/3], Eq. (|i~9"D is accurate to the 
percent level, and this is adequate in comparison with the error of the helioseismological 
determination of c&, see Eq. (|T5|) . 

On the other hand, the density pb at the bottom of the convective zone is an independent 
quantity. Actually, from the above equations p(x) in the convective zone is determined up 
to a scaling factor. The helioseismological determination of pb fixes such a factor. Thus pb 
is also an interesting quantity, which however received little attention in the past. 

We derived pb from U(x) by integrating the hydrostatic equilibrium equation: 



I dp _ 1 

p dx U 

together with 

dM T 



dU GM X 



(20) 



dx 



47rpR Q x 2 . (21) 



Thus pb is sensitive to all the uncertainties of the inversion method, including those of small 
x region. 

Our analysis, see again Table ||, yields: 

R b /R Q = 0.708 - 0.714 (22) 



essentially confirming the result of ||13|| . Note that again the main sensitivity is to the 
starting model. Concerning the density, we find: 

p b = (0.185 - 0.199) g/cm 3 (23) 

Also in this case, a large fraction of the error comes from the starting model. 



V. THE INTERMEDIATE REGION 

The essential output of helioseismology is the reconstruction of the sound speed profile. 
Our discussion is in terms of the related quantity U; as well known, c 2 =TU and below the 
convective zone T = 5/3 with an accuracy of 10~ 3 or better. 

By using the RLS method, it is possible to derive directly the profile of U as a function 
of the radial coordinate throughout all the sun, except for the inner region (x < 0.1). The 
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results (values of U, partial uncertaintes and global errors), are shown in Table [IV| and Fig. 
H] and are summarized in Table [Tl[ 

It is convenient to consider an intermediate solar region: 0.2 < x < 0.65. The upper 
limit is established by requiring that it is well below the transition to the convective zone, 
which we discussed above. The lower limit is chosen so as to exclude the region of energy 
production (see next section). For this region the following comments are relevant: 

a) Each of the individual uncertainties nowhere exceeds 2 °/ 00 - 

b) Uncertainties from the accuracy on the measured frequencies are of minor relevance 
with respect to the residual model dependence and to the sensitivity to the inversion pa- 
rameters. 

c) All in all, even with the most conservative estimates, the helioseismological determi- 
nation is extremely accurate: \AU/U\ < 5°/ 00 throughout the explored region. 



VI. THE ENERGY PRODUCTION REGION 

As well known, most of the energy and of solar neutrinos originate from the innermost 
part of the sun. According to SSM calculations, see e.g. Refs. |3|j3^], about 94% of the solar 



luminosity and 93% of the pp neutrinos are produced within x < 0.2, the region which we 
analyse in this section. 

Our results, summarized in Fig. H and Table f|, deserve the following comments: 

i) In the region 0.1< x <0.2, where the RLS method holds, the error intrinsic to the inversion 
method is dominant, and the global uncertainty, which worsens as one is going deeper into 
the sun, reaches ± 1%. 

ii) At even smaller radii, where the hybrid method is used to determine the function U(x), 
a large fraction of the error comes from the choice of Xf and of xq, whereas the residual 
dependence on A and np is rather weak. The result is also essentially unsensitive to the 
starting model. On the other hand, errors on frequencies become important. This obviously 
corresponds to the fact that p-modes do not penetrate in the solar core , and consequently 
the information one can extract from available experimental results is limited. 

hi) We remind that the production of 7 Be neutrinos is limited to a small region close to the 
solar center, the maximal production occurring at: 

x Be = 0.06 (24) 

according to ||. The 8 B neutrinos are generated even closer to the center, the maximal 
production being at M: 

x B = 0.04. (25) 

(The exact locations of these maxima are somehow model dependent: for instance, in the 
solar model with helium and metal diffusion of Ref. ||34|| , hearafter referred to as BP95, one 
has XBe — 0.055 and xb — 0.045). One sees from Table || that U(x Be ) and U(x B ), as derived 
from the hybrid method, are globally known with an accuracy of about 2%. The average 
values around xs e and (xb) given by the SOLA method are in agreement, within the quoted 
errors, with the results of the hybrid method. 
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VII. HELIOSEISMOLOGY AND STANDARD SOLAR MODELS 



A few recent solar models will now be compared with the helioseismological information 
derived in the previous sections. 

In this respect, the predicted vs. observed properties of the convective envelope are 



particularly interesting, see Table |TXT] and Fig. |3|. All solar models neglecting diffusion 
cannot account for Y p h and for the properties of the bottom of the convective zone (Rb 
and pb). Models with diffusion of helium and heavier elements are much closer to the 
helioseismological determinations. Among these, BP95 and JCD yield values of Y p h, and 
Pb in good agreement with helioseismological predictions. 

For understanding what is going on, let us concentrate on our own solar calculations. 



FRANEC96 [35| is our most recent solar model, including diffusion of He and heavier el- 



ements and using the Livermore EOS and the new Livermore opacities for 19 elements 



according to the composition of | 3~Tfl . If diffusion is switched off, the resulting model 
(FRANEC96-ND) is grossly inconsistent with helioseismology. FRANEC96 looks much bet- 
ter in this respect. Hovewer it slightly underestimates the depth of the convective zone; as 
this is too shallow, the density at the bottom is too small (and correspondingly the sound 
speed is somehow underestimated). If we use the older Livermore opacities calculated for 12 
elements, the resulting model (FRANEC96-OLD) is in good agreement with helioseismology. 

In Fig. [| we show the difference between U as predicted by selected solar models and 
the helioseismological determination, normalized to this latter. We remark that BP95 and 
JCD model are in agreement with helioseismology everywhere. 



SUN24, the model without diffusion of Ref. |^3| yields a good profile of U(x). We 
remind, hovewer, that it fails in predicting the properties of the convective envelope. This 
shows that the two approaches (profile of U and properties of the convective envelope) are 
complementary and both important. 

FRANEC96 underestimates U by about 1% near the transition between the convective 
and radiative region, as already discussed. On the other hand, FRANEC96-OLD, as well as 
other models using older versions of the Livermore opacities (BP95, JCD and SUN24) look 
better. This does not imply some deficiency of the latest opacities. The important point, 
however, is that helioseismology can discriminate between models based on opacities which 
differ just by a few percent. This is best seen in Fig.[| where the profile of U according 
to FRANEC96 and FRANEC96-OLD are compared, together with the relative differences 
Ak/k between opacities, these latter being calculated at any point x for the same value of 
p, T and chemical composition (given by FRANEC96-OLD). 

In conclusion, it seems to us that all these models with diffusion essentially pass the 
helioseimological tests. Similar conclusions have been recently reached in Ref. |48[], where 
the inconsistency of mixed solar models with helioseismology is also discussed. 

VIII. CONCLUDING REMARKS 

We summarize here the main points of this paper: 

• An extensive an possibly exhaustive investigation of uncertainties in helioseismic de- 
terminations of solar properties has been performed. 
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• The accuracy on quantities characterizing the convective envelope (Y p h, Rb and p&) 
has been studied, see Table [□], as well as that on the squared isothermal sound speed, 
U = P/p, along the solar profile, see Fig. |2|. 

• Only recent standard solar models including diffusion of He and heavier elements 
predict properties of the convective envelope in agreement with helioseismology, see 
Table [TTT| and Fig. g. 

• The same models also predict for U(x) values consistent with the helioseismic deter- 
mination all over the solar profile, see Fig. |j. 

Solar models built ad hoc so as to solve the solar neutrino problem should be at least as 
successful. 
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TABLES 



TABLE I. Helioseismological determinations of the photospheric He abundance, Y p h- 



Ref. 




Dappen (1991) @ 


0.268 ± 0.010 


Vorontsov fl99L) HTl! 


250 + 010 


Dziembowski (1991) @ 


0.234 ± 0.005 


Kosovichev (1992) fTl] 


0.232 ± 0.006 


Cox (1993) || 


0.240 ± 0.005 


Dziembowski (1994) p8[ 


0.24295 ± 0.0005 


Hernandez (1994) || 


0.242 ± 0.003 


Antia (1994) @ 


0.252 ± 0.003 


Gough (1996) gl| 


0.248 ± 0.005 


RVCD96 ffl 


0.250 ± 0.005 



TABLE II. For the indicated quantities Q we present the helioseismological reference values 
Qo, the relative partial uncertainties (SQ/Q)k corresponding to experimental errors (exp), model 
dependence (mod), and to the parameters used in the inversion method (A, rip, Xfa and x a ). Global 
errors AQ/Q, estimated according to the statistical (sta) and the conservative (con) approach, see 
Eqs. (|l~0|) and (|TTj) , are also shown. All uncertainties and errors are in °/ 00 . In the fifth and sixth 
row, for the quantity U = Pj p the values of the partial uncertainties and of the global errors are 
the maxima in the indicated interval. In the last two rows the results on U at points representative 
of the 7 Be and 8 B neutrino production, Eqs. (^) and (|25|), have been derived by using the hybrid 
method. 



Q 



exp 



(f) 



mod 



(SQ\ [SQ\ (SQ\ [SQ\ 
\Q )\ \ Q)n F \ QJx fit \ QJx 



(AQ\ 
V Q Jsta 



V Q A 



Yph 

Rb/Re 
Pb [g/cm 3 ] 
U(<d.2 <x< 0.65 ) 
1/(0.1 < x < 0.2 ) 
U(x Be ) [10 15 cm 2 s~ 2 ] 
U(x B ) [10 15 cm 2 s" 2 ] 



0.249 
0.711 
0.192 



1.56 
1.55 



2.4 
0.1 
3 
< 1 
<1 
2.7 
4.8 



27 

4 

16 
~Y~ 

1.6 
To" 

1.3 



6 

0.03 
4 

1.3 

2.5 
1.7 

<1 



7 
0.01 
6.4 
1.4 
3.0 
2.1 
1.5 



2.5 



1.5 
3.7 



6.6 
11 



14 

2 

9.4 

TT 

2.3 
T3~ 
6.6 



42 
4 
37 

"T" 

9.4 
"l8~ 

24 
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TABLE III. Helioseismological determinations and solar model predictions for the properties 
of the convective envelope 



model 



Ref. 



Rb/Rd 



© 



Models without diffusion: 

FRANEC96-ND 
BP95-ND 
CGK89-ND 
P94-ND 
RVCD96-ND 
DS96-ND 
BCDSTT-ND 
TCL 
SUN24 
SUN241 



m 



m 
m 



m 



m 



9 



0.728 
0.726 
0.714 
0.726 
0.725 
0.731 
0.721 
0.725 
0.716 * 
0.714 * 



Models with He diffusion: 

BP92 ! I 

P94 

BCDSTT 



4| 



0.707 
0.710 
0.707 



Y 



ph 



0.261 
0.268 
0.291 
0.270 
0.278 
0.285 
0.276 
0.271 
0.283 * 
0.282 * 



0.247 
0.246 
0.246 



p b [g/cm 3 



0.156 
0.157 



0.166 

0.171 * 
0.166 * 
0.184 * 
0.186 * 



0.197 



0.199 



Models with He and Z diffusion: 

FRANEC96 || 

FRANEC96-OLD 

BP95 g 

CGK89 g| 

P94 g 

RVCD96 H 

DS96 || 

JCD m 



0.716 
0.713 
0.712 
0.721 
0.712 
0.716 
0.713 
0.711 



0.238 
0.242 
0.247 
0.256 
0.251 
0.258 
0.231 
0.245 



0.181 
0.187 
0.187 



0.188 
0.190 



Helioseismology 



0.708-0.714 



0.238-0.259 



0.185-0.199 



private comunication by the authors 
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TABLE IV. We present as a function of x = R/Rq the helioseismological (best) values for 
U = P/p (Uq), the relative partial uncertainties (8U/U)k corresponding to experimental errors 
(exp), model dependence (mod), and to the parameters used in the inversion method (A, np, xju 
and x Q ). Global errors AU/U, estimated according to the statistical and the conservative approach, 
sec Eqs. © and (0), are also shown. 



X 


U Q [c.g.S.] 


\ u J eX p 


(EL) 

V U J mod 




w 


A 




( 5U\ 
\ — ) r 


F 






V U J sta 


\ u J con 


0.0084 


1 


539E+15 


6 


829E- 


03 


1 


170E- 


03 


1 


589E- 


04 


1 


004E- 


03 


1.315E-02 


1.591E-02 


1 


090E- 


02 


3 


823E- 


02 


0.0099 


1 


540E+15 


6 


792E- 


03 


1 


169E- 


03 


1 


617E- 


04 


1 


013E- 


03 


1.296E-02 


1.583E-02 


1 


080E- 


02 


3 


792E- 


02 


0.0116 


1 


540E+15 


6 


740E- 


03 


1 


169E- 


03 


1 


733E- 


04 


1 


027E- 


03 


1.269E-02 


1.571E-02 


1 


067E- 


02 


3 


751E- 


02 


0.0136 


1 


542E+15 


6 


670E- 


03 


1 


233E- 


03 


1 


901E- 


04 


1 


046E- 


03 


1.233E-02 


1.555E-02 


1 


050E- 


02 


3 


702E- 


02 


0.0160 


1 


543E+15 


6 


571E- 


03 


1 


242E- 


03 


2 


126E- 


04 


1 


071E- 


03 


1.183E-02 


1.533E-02 


1 


026E- 


02 


3 


626E- 


02 


0.0188 


1 


544E+15 


6 


438E- 


03 


1 


424E- 


03 


2 


434E- 


04 


1 


106E- 


03 


1.116E-02 


1.503E-02 


9 


940E- 


03 


3 


540E- 


02 


0.0221 


1 


546E+15 


6 


257E- 


03 


1 


545E- 


03 


2 


855E- 


04 


1 


154E- 


03 


1.025E-02 


1.462E-02 


9 


511E- 


03 


3 


411E- 


02 


0.0260 


1 


549E+15 


6 


010E- 


03 


1 


638E- 


03 


3 


426E- 


04 


1 


219E- 


03 


9.045E-03 


1.407E-02 


8 


945E- 


03 


3 


232E- 


02 


0.0305 


1 


551E+15 


5 


681E- 


03 


1 


590E- 


03 


4 


825E- 


04 


1 


305E- 


03 


7.473E-03 


1.332E-02 


8 


216E- 


03 


2 


985E- 


02 


0.0357 


1 


554E+15 


5 


249E- 


03 


1 


449E- 


03 


6 


643E- 


04 


1 


419E- 


03 


5.490E-03 


1.234E-02 


7 


324E- 


03 


2 


661E- 


02 


0.0416 


1 


557E+15 


4 


698E- 


03 


1 


348E- 


03 


8 


941E- 


04 


1 


566E- 


03 


3.107E-03 


1.109E-02 


6 


322E- 


03 


2 


271E- 


02 


0.0482 


1 


560E+15 


4 


026E- 


03 


1 


346E- 


03 


1 


172E- 


03 


1 


746E- 


03 


5.675E-04 


9.564E-03 


5 


344E- 


03 


1 


842E- 


02 


0.0555 


1 


562E+15 


3 


248E- 


03 


1 


408E- 


03 


1 


488E- 


03 


1 


957E- 


03 


2.164E-03 


7.783E-03 


4 


578E- 


03 


1 


805E- 


02 


0.0634 


1 


563E+15 


2 


407E- 


03 


1 


472E- 


03 


1 


823E- 


03 


2 


189E- 


03 


4.245E-03 


5.833E-03 


4 


127E- 


03 


1 


797E- 


02 


0.0717 


1 


561E+15 


1 


579E- 


03 


1 


473E- 


03 


2 


144E- 


03 


2 


429E- 


03 


5.052E-03 


3.853E-03 


3 


726E- 


03 


1 


653E- 


02 


0.0805 


1 


557E+15 


9 


057E- 


04 


1 


542E- 


03 


2 


406E- 


03 


2 


661E- 


03 


4.142E-03 


2.035E-03 


3 


056E- 


03 


1 


369E- 


02 


0.0896 


1 


549E+15 


6 


468E- 


04 


1 


549E- 


03 


2 


550E- 


03 


2 


852E- 


03 


2.977E-03 


6.362E-04 


2 


584E- 


03 


1 


121E- 


02 


0.0992 


1 


538E+15 


7 


289E- 


04 


1 


496E- 


03 


2 


495E- 


03 


2 


962E- 


03 


1.809E-03 


3.901E-06 


2 


293E- 


03 


9 


494E- 


03 


0.1093 


1 


523E+15 


8 


151E- 


04 


1 


576E- 


03 


2 


207E- 


03 


2 


993E- 


03 


7.635E-04 


0.000E+00 


2 


095E- 


03 


8 


354E- 


03 


0.1200 


1 


501E+15 


9 


119E- 


04 


1 


553E- 


03 


1 


755E- 


03 


2 


977E- 


03 


8.726E-05 


0.000E+00 


1 


949E- 


03 


7 


284E- 


03 


0.1315 


1 


473E+15 


9 


455E- 


04 


1 


441E- 


03 


1 


178E- 


03 


2 


886E- 


03 


0.000E+00 


0.000E+00 


1 


781E- 


03 


6 


452E- 


03 


0.1440 


1 


438E+15 


8 


597E- 


04 


1 


313E- 


03 


5 


731E- 


04 


2 


690E- 


03 


0.000E+00 


0.000E+00 


1 


583E- 


03 


5 


435E- 


03 


0.1575 


1 


394E+15 
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FIGURES 



FIG. 1. Helioseismological determinations of the photospheric Helium abundance, Y p /j. Same 
notation as in Table III. The allowed region according to the present work is between the dashed 
lines. 



FIG. 2. The estimated global relative uncertainty on U = P/p, according to the conservative 
approach, Eq. ( |Tl| ) (thick line), and to the statistical approach, Eq. ( |i~0| ) (thin line). 

FIG. 3. Helioseismological determinations and solar model predictions of properties of the 
outer sun. The box defines the region allowed by helioseismology. Open circles denote models 
without diffusion, squares models with He diffusion, full circles models with He and heavier elements 
diffusion. 



FIG. 4. The difference between U as predicted by selected solar models, U mo d and the helio- 
seismological determination, U@, normalized to this latter. The allowed region is that within the 
thick lines, corresponding to (4r) ■ SUN24 is the "model 0" of Ref. pi; FRANEC96 is the 

V u I con 

"best" model with He and heavier elements diffusion of Ref. [35|; BP95 is the model with metal 
and He diffusion of Ref. @; JCD is the "model S" of Ref. pi. 



FIG. 5. a) The profiles of U according to FRANEC96 and FRANEC96-OLD (full and dashed 
lines, respectively), b) The relative differences Ak/k between 19- and 12- elements opacities, along 
the solar profile. 
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This figure "fig2.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9612053vl 



This figure "fig4.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9612053vl 
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